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We study chemotaxis of a single E.coli bacterium in a medium where the nutrient chemical is 
also undergoing diffusion and its concentration has the form of a Gaussian whose width increases 
with time. We measure the average first passage time of the bacterium at a region of high nutrient 
concentration. In the limit of very slow nutrient diffusion, the bacterium effectively experiences a 
Gaussian concentration profile with a fixed width. In this case we find that there exists an optimum 
width of the Gaussian when the average first passage time is minimum, i.e., the search process is 
most efficient. We verify the existence of the optimum width for the deterministic initial position 
of the bacterium and also for the stochastic initial position, drawn from uniform and steady state 
distributions. Our numerical simulation in a model of a non-Markovian random walker agrees well 
with our analytical calculations in a related coarse-grained model. We also present our simulation 
results for the case when the nutrient diffusion and bacterial motion occur over comparable time- 
scales and the bacterium senses a time-varying concentration field. 

PACS numbers: 05.40.Jc, 05.10.Gg, 87.17.Jj 


I. INTRODUCTION 

In a wide variety of physical systems, search process plays an important role [lj. Examples can be found in systems 
such as animals searching for food 0, proteins searching for the binding site on DNA Q, or, diffusion-limited reactions 
Q. A search process is often characterized by the first passage time, which is defined as the time taken to reach 
the target for the first time (or complete the search process). An efficient search process corresponds to a short first 
passage time. Therefore it is crucial to determine how the first passage time depends on the system parameters. 
The most efficient search strategy is often determined by looking into the minimum of the first passage time in this 
parameter space d, 0. 

We consider the search process in one of the most well studied biological systems, viz ., E.coli chemotaxis, which 
describes the motion of E.coli bacteria in response to a chemical concentration gradient [ 3 . When such bacteria are 
placed in an inhomogeneous concentration of a nutrient, they show a tendency to migrate towards the nutrient-rich 
region [&, Q. We ask the questions: How long does it take for the bacteria to find the most favorable region, and 
under what conditions is this search process most efficient? 

The motion of E.coli takes place in two different modes: run and tumble. During a run the bacteria move along a 
given direction with a fixed velocity v ~ 20 pan/s, and during a tumble they do not undergo appreciable displacement 
but change their orientation randomly. In a homogeneous medium, the average run duration is about Is, and at 
the end of one run the bacteria go into the tumbling mode, which lasts for about 0.1s before another run in a new 
direction starts [ToL fllj . In the presence of an inhomogeneous nutrient concentration in the medium, the small size 
(~ 2 fim) of an E.coli cell prevents it from directly sensing the spatial gradient of the concentration field. Therefore, 
to navigate to the nutrient-rich region, the bacteria rely on temporal integration and modulate their run durations 
in different directions in the medium via a memory kernel, shown in Fig. □ El- The kernel effectively compares the 
concentration experienced in the recent past to that in distant past, and if the difference is positive (negative) the 
run duration in the current direction is extended (shortened). Following the above description, the motion of a single 
bacterium in a spatially varying chemical environment is modeled as a non-Markovian random walker with run and 
tumble modes, and the switching rates between these modes depend on the nutrient concentration along its recent 
trajectory HU 17 g| . 

In many physical situations, the chemotaxis motion may take place in a medium where the nutrient is also undergoing 
diffusion. Imagine a situation where a puff of nutrient is injected into the medium such that immediately after 
injection, all the chemical is concentrated in a narrow region. As time goes on, this chemical spreads over the medium 
via diffusion and at any given time its concentration profile has the form of a Gaussian. At a certain stage of this 
nutrient diffusion, when the nutrient has already spread through some distance in the medium, a bacterium is released 
somewhere in the medium which would perform chemotaxis. If the time scale of nutrient diffusion is much longer 
than that of bacterial motion, the bacterium would effectively experience a Gaussian concentration profile of fixed 
width. When the nutrient diffusion and bacterial motion occur over comparable time scales, then the concentration 
sensed by the bacterium will be time dependent and can be described as a Gaussian whose width keeps increasing 
(and peak height keeps decreasing) with time. The region around the peak of the Gaussian profile, where the nutrient 
concentration is highest, is therefore the most favorable region for the bacterium. 
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FIG. 1: (Color online) Bilobe response function of wild-type E. coli used in our simulations. For computational simplicity, we 
have used a discrete sampling of the experimental data presented in [2l| instead of working with the complete data set. This 
did not affect our conclusions. 



In this paper, we study chemotaxis motion of a single bacterium in the presence of a Gaussian concentration field of 
the nutrient, using the model of a non-Markovian random walker described above. As considered in the experiments 
involving bacteria in a micro-fluidic channel or capillary assay HM1 , we consider the motion of the bacterium in 
one dimension. We measure the average first passage time of the bacterium at the neighborhood of the Gaussian 
peak using Monte Carlo simulation. In the limit of slow nutrient diffusion, if the chemotaxis starts at a stage when 
the width of the Gaussian profile is er, then this width does not change appreciably during the search process. The 
average first passage time in this case shows a minimum as a function of er. In other words, there is an optimum 
value of er for which the search process is most efficient. This finding is interesting since a is a parameter that 
can be easily tuned in an experiment, and our study shows that when er is set at a special value, the bacterium 
becomes the most efficient searcher and is able to find its favorable region in the shortest possible time. We also 
perform analytical calculations of mean first passa ge t ime within a coarse-grained model which allows an approximate 
Markovian description of the bacterial motion [lj| [iji] . We find reasonably good agreement between our analytics and 
numerics. We consider a deterministic initial condition as well as a stochastic initial condition, drawn from uniform 
and steady state distribution. In all the cases our numerical simulation and analytical calculations show the existence 
of an optimum width of the nutrient concentration profile when the search process is most efficient. 

In the case when the time scale of nutrient diffusion is comparable to that of bacterial motion, the bacterium 
experiences a time-varying concentration profile, a Gaussian whose width increases with time. The search process 
now crucially depends on the nutrient diffusivity, as well as the extent of spread of the nutrient in the medium at 
the onset of chemotaxis. Our simulations show that if the chemotaxis starts at an early stage of nutrient diffusion, 
when the width of the Gaussian is still small, then the mean first passage time shows a minimum as a function of the 
nutrient diffusivity. But if the chemotaxis starts at a late stage, when the nutrient has already spread considerably in 
the medium and the width of the Gaussian profile is large, the mean first passage time increases monotonically with 
nutrient diffusivity. 

In the next section, we present our model and summarize earlier results. Sections 3 and 4 contain our analytical 
and numerical results, respectively, for a Gaussian nutrient concentration. In section 5 we present our results for a 
time-varying concentration field. The conclusion is presented in section 6. 


II. DESCRIPTION OF THE MODEL AND EARLIER RESULTS 
A. Simulation in a model of a non-Markovian random walker 

Following [THL ITii [20j , we model the motion of a single bacterium in one dimension as a non-Markovian random 
walker whose dynamics is governed by runs and tumbles. During a run the bacterium moves along one particular 
direction with a fixed velocity. The duration of a run is a stochastic variable and follows a Poissonian distribution 
with a mean of Is in a homogeneous medium. At the end of a run, the bacterium goes into a tumbling mode, in 
which it rotates about itself in a random fashion, without much net displacement, before it starts running again in 
a new direction. Because the average tumble duration is much smaller than that of runs, tumbles are modeled as 
instantaneous events which allow the bacterium to change its direction between two successive runs. In our present 
one-dimensional model, the probability that the run direction is changed (reversed) is denoted as q. In the presence 
of a nutrient concentration gradient in the medium, the tumbling rate depends on the recent history. The probability 
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that a running bacterium tumbles during a time interval [t, t + dt} is then given by 

J dt'R(t - t')c(O) , (1) 

where r is the mean run duration in a homogeneous environment, c(t') is the concentration experienced at a past 
instant t' < t, and R(t) is the response kernel. R(t) contains information about the signaling pathway present inside 
the bacterial cell, and it was measured experimentally for wild-type bacteria in |12l |2l| . R(t) was shown to have a 
bilobe shape, with a relatively sharp positive lobe at smaller t and a shallow negative lobe at larger t that vanishes 
for t > 4s (see Fig. [[]). The area under the positive lobe is roughly equal to the area under the negative lobe, and this 
ensures that the response kernel merely measures the concentration gradient and is insensitive to any overall change 
in the background concentration level. Because of this property, it is called an adaptive response kernel. 

We are interested in the linear response regime, where the integral in Eq. Q] is much less than unity. Within this 
linear theory, we can decompose the above bilobe response kernel into (5-function response kernels of suitably chosen 
amplitudes, and from a superposition of the solutions for these (5-function kernels, we can obtain the solution for the 
full bilobe response. Hence we first consider R(t) = aS(t — A) and analyze this case in detail, where we keep terms 
only up to first order in a. Later, we generalize our results for the full response kernel. 


B. Analytical calculation in a coarse-grained Markovian description 


In an earlier study fl5| a simple coarse-grained model was proposed for describing a single bacterium in a concen¬ 
tration gradient in one dimension. The bacterium was assumed to be confined in a one-dimensional box of length 
L with reflecting boundary walls. Although the dependence of the tumbling rate on the past trajectory makes the 
underlying run-tumble motion non-Markovian, one can still expect that at a coarse-grained level, a Markovian de¬ 
scription might be possible. For this purpose, one can coarse-grain over a time scale which is much larger than the 
typical run duration r and assume that the average bacterial density within the spatial resolution of coarse-graining 
vt has a Markovian dynamics. For the time evolution of this coarse-grained density P{x, t ) at point x , at time t, the 
following Fokker-Planck equation was formulated in [l5| : 

d t P(x,t) = -d x [V(x)P(x,t) - d x (D(x)P{x,t))\ , (2) 

which is the equation for a random walker with position-dependent drift and diffusion. In [l6| the past memory of the 
cell trajectory was included in the description by introducing additional ‘internal variables’, and the resulting process, 
characterized by a larger number of degrees of freedom, now becomes Markovian. Starting from this Markovian 
description, using the homogenization method jH, [HI the same coarse-grained equation as above was obtained [l6|. 

The chemotactic drift velocity V(x) and the diffusivity D(x) in Eq. [2] depend on the nutrient concentration profile 
c( x), and the dependence can be derived from the microscopic dynamics. In an earlier calculation by de Gennes 
an approximate expression for the drift velocity was obtained within the simplifying assumption that whenever a 
running bacterium tumbles, its past memory is lost. Considering the response function R(t) = aS(t — A), the 
resulting expression 
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V(x)=a——e ~d x c(x ) (3) 

Zq 

was found to show good agreement with the simulation results [151]. To calculate the diffusivity D{x ), the effective 
tumbling frequency was calculated within the coarse-grained model by averaging over a population of non-interacting 
bacteria within the coarse-graining length scale. Although the tumbling frequency of a single bacterium depends 
on the details of its past trajectory, this dependence is lost when averaged over a large number of bacteria with all 
possible run directions. The average tumbling frequency at a position x can be shown to be [1 — q;c(x)]/t, from which 
the diffusivity turns out to be 0 

D(x) = ^-ll+a c (x)}. (4) 

Using Eqs. [3] and HI it can be easily shown from Eq. [2] that in the steady state the bacterial density P{x) has the 
form 


P(x) = P 0 + <xP 0 (e t - 1) 


^c(x) — Pq J c(x)da 


( 5 ) 
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where Pq = l/L and the last term takes care of the normalization. 

In the next section, we use the above coarse-grained model to calculate the mean first pa ssag e time of the bacterium 
at the nutrient-rich region. Note that in contrast to the steady state behavior, studied in [15|], we study first-passage 
properties, away from steady state. 


III. ANALYTICAL CALCULATION OF FIRST PASSAGE TIME 


Let P(y,t |x,0) be the conditional probability to find the bacterium at position y at time t, given that it started at 
x at t = 0. This conditional probability follows the backward Fokker-Planck equation: 

d t P{y,t\x,0) = V(x)d x P(y,t\x,0) + £>(x)<9 2 P(i/, t\x, 0). (6) 

To find out the first passage time at a certain point x$ (which we call the target), we consider an absorbing boundary 
condition at xq, in addition to the reflecting boundary at x = 0 and x = L. Without any loss of generality, we 
perform all our measurements for x < xq. The survival probability G(x , t ) that starting from an initial position x, the 
bacterium will not reach the target until time t can be written as G(x,t) = fg° dyP(y,t\x,0). From Eq. [G]it follows 
that G(x, t) satisfies the following equation: 

d t G{x, t) = V{x)d x G(x, t) + D(x)d 2 x G{x, t), (7) 


with the initial condition G(x, 0) = 1. The reflecting and absorbing boundary conditions are implemented as 
d x G(0,t) = 0 and G(xo,i) = 0, respectively. 

By definition, G(x, t ) is the probability that the first passage time is larger than t, and hence the first passage 
time distribution is given by — d t G(x , t). Mean first passage time T(x) = — / 0 °° dt t d t G(x, t) = / 0 °° dt G(x, t), which 
follows the equation 


V{x)d x T(x) + D(x)d 2 T(x) = — 1. 


The solution of this equation has the form 



( 8 ) 

(9) 


where ip(x) = exp [/ c) l dx'V {x') / D{x')]. Now, using Eqs. [Hand ED and keeping terms only upto first order in a, one 
can write ip{x) = exp [ae -29A / r {c(x) — c(0)}] = 1 + ae~ 2qA t T [c{x) — c(0)]. After a few steps of simple algebra we 
finally have 


T{x)= 


- - a(l - e~ 2qA / T ) f*° dy / 0 y dz c(z) 
—ae~ 2qA / T f*° dy y c{y)\, 


( 10 ) 


which can be written in the form T(x ) = Tq(x) + aTi(x), where T$(x) stands for the mean first passage time for 
an ordinary Brownian motion in the absence of any concentration gradient and T\ (x) gives the first order correction 
term when position-dependent drift and diffusion are present due to a spatially varying concentration field c(x) [24| . 
In the rest of this work we focus on T\ (x). 


For a Gaussian concentration profile c(x) = 


exp[—(x — x) 2 /2cr 2 ] 

V2* 


the drift velocity V(x) and diffusivity D(x) show 


rapid variation close to the peak at x. In our coarse-grained description, which allows for analytical treatment, we 
deal with length scales much larger than the mean free path of the bacteria, and any spatial variation that occurs 
over a smaller length scale must be neglected in our coarse-grained model. When a is not too large, the variation 
of V(x) and D(x) around the peak is too fast to be considered in our coarse-grained formalism. Because of this, we 
choose the target position xq slightly away from the peak such that both the initial position x and the target lie on 
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the same side of the peak. For our choice of x < Xq < x we use the Gaussian c(x) in Eq. [lllland get 
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(ii) 


So far we have considered the first passage time with a deterministic initial condition, where the bacterium always 
starts from a fixed initial position x. Now we consider the case of stochastic initial position when x can take any 
value within the interval 0 < x < Xq; that is, the initial position can lie anywhere between the left boundary and the 
target, with a certain distribution function. We consider two specific cases: (*) when x follows a uniform distribution 
Pq and (ii) when x is drawn from the steady state distribution P(x) in Eq. [5j In the first case, the a order correction 
in the first passage time can be obtained by simply integrating Eq. [IT] over x: 
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For case (ii) the initial position x follows the steady state distribution in Eq. [5] The mean first passage time is then 
written as 
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For a Gaussian form of c(x) the a order term becomes 
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where To (a;) and Ti(a:) are defined in Eqs. [TUI and HT1 After straight forward algebra the a order term in first passage 
time with the steady state initial condition becomes 
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FIG. 2: (Color online) Results for the mean first passage time with a fixed initial position. The main plot shows T\(x) (in 
seconds) as a function of standard deviation a (in (im) of the Gaussian nutrient concentration field with R(t) = ad(t — A) 
and A = 0.5s, x = 300/rm. These data have been averaged over 10 7 different trajectories. The top right and bottom right 
insets show the variation of the optimum width a* (in fim) as a function of A (in seconds) and the initial position x (in 
respectively. The error bar lies within 0.1 fim. The top left inset shows Ti(x) (in seconds) vs er (in fim) variation, averaged over 
10' trajectories, for the bilobe response kernel, shown in Fig. |T| The discrete symbols correspond to simulations and the solid 
lines correspond to analytical calculations. Here L = 1000/xm, xo = 490 fim, x = 500 fim, q = 0.5, r = Is. 


The results in Eqs HU HU and HU are for an impulse response kernel R(t) = aS(t — A), and these can be easily 
generalized for any arbitrary response function. In the next section, we measure the first passage time in simulation 
and compare the result with above analytical calculation. 


IV. SIMULATION RESULTS ON FIRST PASSAGE TIME 

We perform a Monte Carlo simulation on a one dimensional box of size L, with reflecting boundary walls at the 
two ends. In a time interval dt the bacterium moves by a distance vdt. At the end of each time step, we compute the 
tumbling probability, as in Eq. HI For an impulse response function R(t) = aS(t — A), the tumbling probability at 
time t is given by dt/r (1 — ac[x(t — A)]), where x(t — A) is the position of the bacterium at a time A back in the past 
and c[x(t — A)] is the concentration experienced by the bacterium at that past instant of time. At the end of one time 
step the bacterium attempts to tumble with this probability. If the tumbling attempt is unsuccessful, it continues to 
move in the same direction with same velocity v, and if the tumbling attempt is successful, the bacterium changes 
its direction (in this one-dimensional case, it reverses the sign of v ) with probability q. Starting from a given initial 
position x at t = 0, we measure the first passage time at a position Xo and average over different trajectories. In order 
to avoid the region with rapid spatial variation of the concentration field, we consider a target which is one mean free 
path away from the Gaussian peak in the same direction as the starting position of the bacterium: xo = x — vt. 

In Fig. [2] we show the variation of Tj(x) with a (discrete symbols) and compare it with our analytical result in 
Eq. HU (solid lines). We find reasonably good agreement between our simulation and analytical calculation. For very 
small (7 the concentration variation can be perceived only within a very narrow region around the peak x. Hence 
the bacterial trajectory starting from x and ending at xo, which does not cross the peak at x, consists of isotropic 
diffusion for the most part. As a result, in the limit of small a the mean first passage time is given by that for an 
ordinary Brownian motion and is equal to To(x) in Eq. 1101 and the first order term T\(x) goes to zero. Similarly, in 
the limit of very large a the profile is almost flat, and even in this case the motion is close to isotropic diffusion and 
Xj(x) vanishes. Our simulation and analytical calculation are consistent with this simple argument. 

For intermediate er values, Tj(x) must show a non-monotonic variation, since it vanishes for small and large a. 
We find a minimum for Tj(x) at a particular width a*. In other words, there exists an optimal width a* when the 
first passage time at the nutrient-rich region becomes shortest and the bacterium becomes the most efficient searcher. 
The bottom right and top right insets in Fig. [U show the variation of the optimal width a* as a function of the 
initial position x and the memory A of the bacterium, respectively. Note that even in the Markovian limit, when the 
bacterium does not have any memory, and does not accumulate in the nutrient-rich region in the long time limit [l5j , 
its first passage properties still show the existence of an optimal width when the search is most efficient. 
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FIG. 3: (Color online) Mean first passage time with stochastic initial positions. The top row shows data for impulse response 
R(t) = a5(t — A), and the bottom row shows data for the bilobe response. (A) and (B) show data for a uniform initial condition, 
when x can take any value in the range [0, io] with uniform probability Po- (C) and (D) show data for a steady state initial 
condition, when the value of x in the range [0, *o] is drawn from the steady state distribution P{x) in Eq. [o] Here the first 
passage time is measured in the units of seconds, and a is in pan. The other simulation parameters are the same as in the main 
plot of Fig. O These data have been averaged over 10 7 trajectories. The discrete symbols are for numerical data, and the solid 
lines are for analytical calculations. 


For wild-type bacteria, the response kernel has a bilobe shape, and we can reconstruct the kernel as a linear 
superposition of impulse response functions with suitable amplitudes, as shown in Fig. [0 We use this response kernel 
to calculate the mean first passage time for various a. Even for this adaptive bilobe kernel, we find there exists an 
optimal width when the mean first passage time hits a minimum. Our analytical calculations show similar results. 
Interestingly, the value of the optimum width cr* does not change even when the initial position and the target position 
are varied (data not shown here). In other words, a wild-type E.coli bacterium becomes the most efficient searcher 
when placed in an environment with a Gaussian concentration profile of a nutrient with a width of ~ 50pm. 

Instead of starting from a fixed position, when the initial position is a random variable which can choose any value 
between the left boundary wall at x = 0 and the target at x = Xq with a certain distribution, our results show the 
existence of an optimum a that minimizes the first passage time. We have considered initial positions chosen from 
uniform distribution as well as from steady state distribution in Eq. [5j We present our data in Fig. [3j Note that for 
the choice of a steady state initial condition, P) ' does not vanish in the limit of small a but approaches a constant 
value. In other words, even when the width of the nutrient concentration profile is vanishingly small, the first passage 
time of the bacterium is not the same as in a homogeneous medium. In fact, when the system is in the steady state, 
the bacterium has explored the full system and has already experienced the narrow concentration profile present in 
the middle of the box. The steady state measure P(x) is therefore not the same as Po but contains information about 
the narrow concentration field. This gives rise to a non-vanishing a order correction term in the limit a — > 0. 


V. TIME-VARYING CONCENTRATION FIELD 

In this section, we consider the case when the nutrient diffusion in the medium occurs over a time scale comparable 
to that of bacterial motion. The bacterium will then experience a time-varying concentration field. Our analytical 
formalism in section 3 does not work in this case, and we study the system using numerical simulations. The simplest 
description of the nutrient concentration profile can be given by a Gaussian whose width is increasing with time: 

4M) = «p(-2^ 

starts, and T> is the nutrient diffusivity. 

The bacterial motion will depend on ao and D, depending on the time scale t c ~ Uq/T). For f < f c the motion 
depends on a o, and for t t c the motion is mainly controlled by T>. In the limit of very small T>, therefore, one 
would expect the first passage time to be a function of ao alone. In fact this is the limit when the nutrient diffusion 
is very slow, and during the time interval of the first passage at the target, the width of c(x,t) changes very little. 
In this limit, therefore, one expects results similar to those in a static concentration profile. Our simulation data 
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FIG. 4: (Color online) First passage time for time-dependent concentration of the nutrient. (A) shows the variation of T\{x) 
(in seconds) as a function of ao (in /rm) with T> held fixed at 0.01 fj,m 2 /s (red squares), 10/im 3 /s (green circles), and 37 )j,m 2 /s 
(black diamonds). (B) shows T\(x ) (in seconds) vs T) (in /im 2 /s) plot for ao = 30/rm (red squares), 50 fim (green circles) and 
120 fim (black diamonds). The other simulation parameters are the same as in Fig. 0main plot. These data are averaged over 
10 8 trajectories. 


in Fig. |4j4 indeed show that for small T> there is an optimum width ao where T\(x) becomes minimum. As T> 
increases, t c becomes smaller when T\(x) does not show much variation with op and the minimum becomes less and 
less pronounced. In Fig. [4]4 we verify this. 

In Fig. 0J3 we show the variation of Tj(x) against V for fixed op values. For very large T> the Gaussian profile 
quickly flattens out, and the bacterial motion becomes an isotropic diffusion. In this limit T\(x) becomes zero. For 
very small T> values, the limit for a static Gaussian profile is recovered, and (as shown in our data in Fig. [2] main 
plot) T\{x) has a negative value that depends on op. Therefore, for a given op, as V is varied, T\(x) starts from a 
negative value at small T> and becomes zero at large T>. Whether this variation is monotonic or not depends on the 
choice of op. Om data in Fig. 0J3 show that for large ao the variation is monotonic but for small ao a minimum is 
reached at a particular T >; that is, there is an optimum diffusivity of the nutrient when the search is most efficient. For 
our various choice of <to values over a wide range (full data set not presented here), we also notice that an optimum 
diffusivity is observed whenever ao is fixed at a value smaller than a *, the optimum width for the static concentration 
profile (see Fig. [2] main plot). For ao > a*, on the other hand, T\[x) increases monotonically with V. 

The above observation tentatively indicates that it may be possible to describe the results for the time-dependent 
nutrient concentration in terms of a static concentration profile with an ‘effective width’ a e . For a given value of ao 
and T> the width of c(x,t) keeps increasing during bacterial motion: at the start of the motion the width is ao, and 
at the end of the first passage the average width is \J a^ + 2T>T(x). Let us assume that a e is some measure of the 
average or effective width experienced by the bacterium during this process. Obviously, a e is a function of both ao 
and T>: for a fixed ao, as T> is varied, a e ps ao for very small T>, and as T> becomes very large, so does a e . In the 
course of this variation, if a e crosses a*, then T\(x) shows a minimum and if ao > a* such that a e never reaches a* 
(because op can never fall below ao), then T\{x) shows a monotonic increase with T>. The above picture explains our 
numerical data well. However, we would like to mention that we do not yet have any mathematical expression for a e 
in terms of op and T>. It would be interesting to directly verify the mechanism proposed above. 


VI. CONCLUSION 

In this paper, we have considered the chemotaxis motion of a bacterium in a medium where the nutrient is also 
undergoing diffusion and its concentration profile is given by a Gaussian whose width increases with time. We have 
measured the mean first passage time of the bacterium at the neighborhood of the Gaussian peak. In the limit 
when the nutrient diffusion is slow compared to the bacterial motion, the bacterium experiences an effectively static 
concentration profile, a Gaussian with a fixed width, and in this regime we calculate the mean first passage time 
analytically, within a coarse-grained formalism. We find that the mean first passage time shows a minimum as a 
function of the width of the Gaussian, which means that the search process becomes most efficient at a certain 
optimum width. Our numerical simulation matches the analytical result well. 

For a time-dependent concentration profile, i.e., in the regime in which the nutrient diffusion occurs over a time-scale 
comparable to bacterial motion, we find that the first passage time is a function of nutrient diffusivity T> and the 
width op of the Gaussian at the onset of chemotaxis motion. When V is held fixed at a small value, the mean first 
passage time shows a minimum against variation of ao, as in the static case. But for large T> the minimum becomes 
less pronounced. As a function of T>, the mean first passage time shows a minimum if op is held fixed at a small value. 
But no such minimum is observed when a 0 is set at a large value. 
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The existence of an optimal width of the nutrient concentration profile that makes the search process most efficient 
is an interesting result. Apart from the mean first passage time, we have also examined our conclusion by measuring 
the most probable first passage time. Note that for a general process described by an equation of the form given in 
Eq. m the probability distribution for first passage time shows a long tail which makes the mean first passage time 
much larger than the typical (or most probable) one. Our numerical simulations (data not shown here) show that the 
typical first passage time also becomes minimum at a particular width whose value is very close to the one for which 
the mean showed a minimum. 

In [l(| the chemotactic efficiency was characterized by a quantity called ‘uptake’ (defined as the total amount of 
nutrient encountered by the bacterium upto a certain time). It was shown that under harsh environmental conditions, 
when the nutrient is scarce, the long time uptake can be maximized for a particular shape of the response kernel. 
In contrast, we consider a given response kernel and vary the parameter(s) characterizing the concentration gradient 
of the nutrient in the medium to find the minimum first passage time. In other words, for a given response kernel 
we find the optimum environmental conditions for the fastest search process. It is easy to argue that the long time 
uptake does not show any maximum in our case, when the response kernel is fixed and the environmental conditions 
are varied. The long time uptake decreases monotonically as the concentration gradient is decreased. 

Finally, it would be interesting to verify our conclusions in experiments. Recently, E.coli chemotaxis has been 
studied in a micro-fluidic channel whose width is comparable to the bacterial mean free path [l7l-fl9|. In such a 
setup, the motion of the bacterium can be considered to be effectively one-dimensional. It is possible to generate a 
Gaussian chemical concentration profile using techniques of diffusive microfluidics [25j. The motion of the bacterium 
can be tracked to measure its first passage properties. Our model predicts that for a static Gaussian profile of width 
a ~ 50 \xm. wild-type E.coli have shortest first passage time. However, it can be experimentally challenging to verify 
our results for a time-dependent nutrient concentration profile. Most common chemoattractants such as aspartate 
and serine have diffusivity T> ~ 1000 fiirr/s which is much larger than bacterial diffusivity. As a result, the chemical 
diffuses very quickly in the medium, and initially localized concentration quickly flattens out. Thus the bacterium 
experiences a very weak concentration gradient and the chemotactic correction T± to its first passage time may become 
too small for experimental detection. 
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